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Abstract 

Using data collected at the ^(3770) resonance with the CLEO-c detector at the Cornell e+e" 
storage ring, we present improved measurements of the absolute branching fractions of decays 
to K^e'^Ue, vr^e+i/e, K*^e'^i'e, and p^e'^Ue, and the first observation and absolute branching fraction 
measurement of oje'^Ve. We also report the most precise tests to date of isospin invariance 

in semileptonic and decays. 
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The quark mixing parameters are fundamental constants of the Standard Model (SM) 
of particle physics. They determine the nine weak-current quark coupling elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix The extraction of the quark couplings is 
difficult because quarks are bound inside hadrons by the strong interaction. Semileptonic 
decays are the preferred way to determine the CKM matrix elements as the strong interaction 
binding effects are confined to the hadronic current. They are parameterized by form factors 
that are calculable, for example, by lattice quantum chromodynamics (LQCD) and QCD sum 
rules. Nevertheless, form factor uncertainties dominate the precision with which the CKM 
matrix elements can be determined Q. In charm quark decays, however, couplings Vcs and 
Vcd are tightly constrained by the unitarity of the CKM matrix. Therefore, measurements 
of charm semileptonic decay rates and form factors rigorously test theoretical predictions. 

We report herein measurements with the first CLEO-c data of the absolute branching 
fractions of decays to K^e'^i'e, n^e'^i'e, K*^e^Vf.i and p°e"'"z/e, and the first observation 
and absolute branching fraction measurement of uje^Vf.. (Throughout this Letter 

charge-conjugate modes are implied.) We combine these results with the measurements 
of semileptonic branching fractions reported in Ref. 4], which use the same data and 
analysis technique, and test isospin invariance of the hadronic current in semileptonic decays. 

The data were collected by the CLEO-c detector at the ip{?)77Q) resonance, about 40 MeV 
above the DD pair production threshold. A description of the CLEO-c detector is provided 
in Ref. and references therein. The data sample consists of an integrated luminosity of 
55.8 pb~^ and includes about 0.16 million D^D~ events. 

The technique for these measurements was first applied by the Mark III collaboration 
at SPEAR. Candidate events are selected by reconstructing a D~ , called a tag, in the 
following six hadronic final states: i^^vr", K^ti~tt~ , Kgir'n^, K^tt^tt'tt'^ , Kg'n'~TT~7T~^ , and 
K~K~^7i~ . The absolute branching fractions of D'^ semileptonic decays are then measured 
by their reconstruction in the system recoiling from the tag. Tagging a D~ meson in a 
-0(3770) decay provides a with known four-momentum, allowing a semileptonic decay 
to be reconstructed with no kinematic ambiguity, even though the neutrino is undetected. 

Tagged events are selected based on two variables: AE = Ed — -Ebeam, the difference 
between the energy of the D~ tag candidate (Ed) and the beam energy (-Ebeam), and the 
beam-constrained mass Mbc = \J E^^^^/ — \vd\^ !(?■, where is the measured momentum 
of the D~ candidate. Selection criteria for tracks, 7r° and K^^ candidates for tags are described 
in Ref. If multiple candidates are present in the same tag mode, one candidate per tag 
charge is chosen using IS.E. The yields of the six tag modes are obtained from fits to the 
Mbc distributions. The data sample comprises approximately 32,000 charged tags (Table P). 

After a tag is identified, we search for a positron and a set of hadrons recoiling against the 
tag. (Muons are not used as D semileptonic decays at the '?/'(3770) produce low momentum 
leptons for which the CLEO-c muon identification is not efficient.) Positron candidates, 
selected with criteria described in Ref. j^, are required to have momentum of at least 
200 MeV/c and to satisfy | cos^| < 0.90, where d is the angle between the positron direc- 
tion and the beam axis. The efficiency for positron identification rises from about 50% at 
200 MeV/c to 95% just above 300 MeV/c and is roughly constant thereafter. The rates for 
misidentifying charged pions and kaons as positrons averaged over the momentum range are 
approximately 0.1%. Bremsstrahlung photons are recovered by the procedure described in 
Ref. 1^. 

Hadronic tracks must have momenta above 50 MeV/c and | cos 6*1 < 0.93. Identification of 
hadrons is based on measurements of specific ionization [dEjdx) in the main drift chamber 
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Tag Mode 




Yield 


D K^^TT 




2243 ± 51 
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D- K+TT- 




4734 ± 91 


D- K'^gTT- 




3281 ± 94 


D- K^K 




1302 ± 44 


All Tags 




31922 ±219 



TABLE I: Tag yields of the six D hadronic modes with statistical uncertainties. 

and information from the Ring Imaging Cherenkov Detector (RICH). Pion and kaon candi- 
dates are required to have dE/dx measurements within three standard deviations (3. Oct) of 
the expected value. For tracks with momenta greater than 700 MeV/c, RICH information, 
if available, is combined with dE/dx. The efficiencies (95% or higher) and misidentification 
rates (a few per cent) are determined with charged pion and kaon samples from hadronic D 
decays. 

We form 7r° candidates from pairs of photons, each having an energy of at least 30 MeV, 
and require that the invariant mass of the pair be within 3.0o" (a ~ 6 MeV/c^) of the 
known 7r° mass. A mass constraint is imposed when candidates are used in further 
reconstruction. The candidates are formed from pairs of oppositely-charged and vertex- 
constrained tracks having an invariant mass within 12 MeV/c^ (~ 4.5cr) of the known 
mass. We form a K*'^ (p'') candidate from K~ and tt"'" (tt" and vr"*") candidates and require 
an invariant mass within 100 MeV/c^ (150 MeV/c^) of its mean value. The reconstruction of 
uj — > 7r+7r^7r° candidates is achieved by combining three pions, requiring an invariant mass 
within 20 MeV/c^ of the known mass, and demanding that the charged pions do not satisfy 
interpretation as a K^. 

The tag and the semileptonic decay are then combined, if the event includes no tracks 
other than those of the tag and the semileptonic candidate. Semileptonic decays are identi- 
fied using the variable U = £^miss — c|pmiss|, where ^'miss and Pmiss are the missing energy and 
momentum of the D meson decaying semileptonically. If the decay products of the semilep- 
tonic decay have been correctly identified, U is expected to be zero, since only a neutrino is 
undetected. The resolution in U is improved using constraints described in Ref. jj]. Due to 
the finite resolution of the detector, the distribution in U is approximately Gaussian, cen- 
tered at U = with cr ~ 10 MeV. (The width varies by mode and is larger for modes with 
neutral pions.) To remove multiple candidates in each semileptonic mode, one combination 
is chosen per tag mode, based on the proximity of the invariant masses of the Kg, K*^, p^, 
7r°, or a; candidates to their expected masses. 

The yield for each semileptonic mode is determined from a fit to its U distribution, as 
shown in Fig. ^ with all tag modes combined. In each case the signal is represented by a 
Gaussian and a Crystal Ball function 0] to account for initial and final state radiation (FSR). 
The parameters describing the tails of the signal function are fixed with a GEANT-based 
Monte Carlo (MC) simulation 8]. The background functions are determined by a MC sim- 
ulation that incorporates all available data on D meson decays. The backgrounds are small 
and arise mostly from misreconstructed semileptonic decays with correctly reconstructed 
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FIG. 1: Fits (solid line) to the U distributions in data (dots with error bars) for the five 
semileptonic modes: (a) D+ K°e+i^e, (b) L>+ K*°e+i/e, (c) D+ 7r°e+i^e, (d) D+ p°e+i/e 
and (e) — > ooe^Ve- The background contribution is the dashed line only visible in (c) and (d). 
The arrows show the signal region for — > uje^v^,. 



tags j9|. The background shape parameters are fixed, while the background normalizations 
are allowed to float in all fits to the data. 

The mode —>■ uje^Ve has never previously been observed. There are 8 events consistent 
with — >■ uje^Ue in Fig.H(e). The background in the signal region ([—25; +30] MeV in U) 
is estimated to be 0.4 ± 0.2 events. The probability for the background of 0.6 events to 
fluctuate to 8 or more events is 2.4 x 10~^, which corresponds to significance exceeding 5. Oct. 
Therefore, this is the first observation of ue^Ve- 

The absolute branching fractions in Table ITTl are determined using B = A^'signai/e^^tag, 
where Ai'signai is the number of fully reconstructed D^D~ events obtained by fitting the U 
distribution, A^tag is the number of events with a reconstructed tag, and e is the effective 
efficiency for detecting the semileptonic decay in an event with an identified tag. A MC 
simulation where the relative population of tag yields across tag modes approximates that 
in the data is used to determine the efficiency. 

We have considered the following sources of systematic uncertainty and give our estimates 
of their magnitudes in parentheses. The uncertainties associated with the efficiency for find- 
ing a track (0.7%), vr" (2.0% for D+ ue+u^ and 4.3% for D+ n^e+u^) and K^g (3.0%) are 
estimated using missing mass techniques with the data Details on the uncertainties as- 
sociated with positron identification efficiency (1.0%) are provided in Ref. Uncertainties 
in the charged pion and kaon identification efficiencies (0.3% per pion and 1.3% per kaon) 
are estimated using hadronic D meson decays. The uncertainty in the number of tags (1.1%) 
is estimated by using alternative signal functions in the fits to the Mbc distributions and 
by varying the end point of the background function The uncertainty in modeling the 
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Decay Mode e (%) Yield 

D+ K^e+Ve 57.1 ±0.4 
D+ 7r°e+z^e 45.2 ± 1.0 
D+ K*^e+Ve 34.8 ± 0.3 
D+^p°e+i/e 40.0 ±1.1 
D+ uje+Ue 16.4 ± 0.6 



B (%) B (%) (PDG) 



545 ± 24 
63.0 ± 8.5 
422 ± 21 
27.4 ± 5.7 
7 

/.D_2.7 



8.71 ±0.38 ±0.37 
0.44 ± 0.06 ± 0.03 
5.56 ±0.27 ±0.23 
0.21 ±0.04 ±0.01 

o.i6t[!:g^±o.oi 



6.7 ±0.9 
0.31 ±0.15 

5.5 ±0.7 
0.25 ±0.10 



TABLE II: Signal efficiencies, yields, and branching fractions in this work and a comparison to 
the PDG [13] • The first uncertainty is statistical and the second systematic in the fourth column, 
and statistical or total in the other columns. The efficiencies do not include subsidiary branching 
fractions. The D'^ K*^e^Ve yield is reduced by 2.4% in the calculation of the branching 
fraction (see the text for details). 

background shapes in the fits to the U distributions (0.4% to 3.3% by mode) has contribu- 
tions from the uncertainties in the simulation of the positron and hadron fake rates as well 
as the input branching fractions in the MC simulation. The uncertainty associated with 
the requirement that there be no additional tracks in tagged semileptonic events (0.3%) is 
estimated by comparing fully reconstructed DD events in data and MC. The uncertainty in 
the semileptonic reconstruction efficiencies due to imperfect knowledge of the semileptonic 
form factors is estimated by varying the form factors in the MC simulation within their 
uncertainties (1.0%) for all modes except p^e^z/g and uje^Ve-, for these a con- 

servative uncertainty (3.0%) is taken, as no experimental information on the form factors 
in Cabibbo-suppressed pseudoscalar-to-vector transitions exists. The uncertainty associated 
with the simulation of FSR and bremsstrahlung radiation in the detector material (0.6%) 



is estimated by varying the amount of FSR modeled by the PHOTOS algorithm and 
by repeating the analysis with and without recovery of photons radiated by the positron. 
The uncertainty associated with the simulation of initial state radiation (e"^e~ DD'y) is 
negligible. There is a systematic uncertainty due to finite MC statistics (0.7% to 4.0% by 
mode). 

Non-resonant semileptonic decays K^-K^e^Ve are background to K*^e^Ve- 

There is evidence from the FOCUS experiment for a non-resonant component consistent 
with an S-wave amplitude interfering with K*^e~^h'e fl3l|. Its contribution, estimated 

to be 2.4% in this analysis, is subtracted in the calculation of the branching fraction of 
— > K*^e^h'e Systematic uncertainties associated with the subtraction are due to 

imperfect knowledge of the amplitude and phase of the non-resonant component (1.0%), and 
its effect on the reconstruction efficiency (1.5%) The lineshapes for semileptonic modes 
with wide resonances are simulated using a relativistic Breit-Wigner with a Blatt-Weisskopf 
form factor. A systematic uncertainty associated with the K*^ lineshape (1.2%) is assigned 
by comparing the {K^n^) invariant mass distribution in the data to alternative lineshapes 
and the non- resonant contribution. For p^e^Ve-, there is insufficient data to constrain 

the non-resonant background or the resonance lineshape. The systematic uncertainties from 
these two sources are expected to be much smaller than the current statistical uncertainty 
for this mode, and are neglected. 

These estimates of systematic uncertainty are added in quadrature to obtain the total 
systematic uncertainty (Table HH): 4.2%, 5.6%, 4.1%, 6.2%, and 7.8% for D+ K^e+Ue, 
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— s> Ti^e^Vei K*'^e'^i>e, p^e^Ve-, and — > uje^Ve-, respectively. 

We now discuss the results presented in this Letter and the semileptonic study in 
Ref. The measured equality of the inclusive semileptonic widths of and mesons 
demonstrates that the source of the lifetime difference between them is attributable to 
differences in the hadronic widths. The widths of the isospin conjugate exclusive semileptonic 
decay modes of the and are related by isospin invariance of the hadronic current. 
The results obtained here and in Ref. allow the most precise tests so far. 

The ratio rl^+^^Oet^") is expected to be unity. The world average value is 1.35±0.19 ■ 



Using our results and the lifetimes of the and we obtain: r^^+^^Oe+l^'') = 

1.00 ± 0.05(stat) ± 0.04(syst). The result is consistent with unity and with two recent less 
precise results: a measurement from BES II using the same technique and an indirect 
measurement from FOCUS Ratios of isospin conjugate decay widths for other 

semileptonic decay modes are given in Table ITTTl 



Ratio Measured Value 

T{D^ K-e+iy)/T{D+ K'^e+u) 1.00 ± 0.05 ± 0.04 

r(DO ^ K*-e+iy)/T{D+ K*^e+v) 0.98 ± 0.08 ± 0.04 
r(.D° ^ P'e+v)/[2 ■ T{D+ ^ p^e+u)] 1.2^°;^ ± 0.1 



TABLE III: Ratios of semileptonic decay widths of and mesons. The uncertainties are 
statistical and systematic. In each case the ratio is expected to be unity. 

As the data are consistent with isospin invariance, the precision of each branching fraction 
can be improved by averaging the and results for isospin conjugate pairs. The isospin- 
averaged semileptonic decay widths, with correlations among systematic uncertainties taken 
into account, are given in Table HVl 



Decay Mode 
D e+Ue 

D ^K* e+Ve 



r (10-^ -ps-i) 

8.38 ±0.20 ±0.23 
0.68 ±0.05 ±0.02 
5.32 ±0.21 ±0.20 
0.43 ±0.06 ±0.02 



TABLE IV: Isospin-averaged semileptonic decay widths with statistical and systematic uncer- 
tainties. For Cabibbo-suppressed modes, the isospin average is calculated for the using 
r(DO) = 2-r(L>+). 



The ratio of decay widths for D — >■ ire^u and D — > Ke'^u provides a test of the LQCD 
charm semileptonic rate ratio prediction Using the isospin-averaged results in Table HV] 
we find ^ff^^^'^f^)^ = (8.1 ± 0.7(stat) ± 0.2(syst)) x lO^^, consistent with LQCD and two 

recent results ji^, 2l||. Furthermore, the ratio ^-^(^j^l^^+u) predicted to be in the range 
0.5 to 1.1 (for a compilation see Ref. Using the isospin averages in Table HV| we find 

^rloZ^Ke+u) = 0-63 ± 0.03(stat) ± 0.02(syst). 
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Finally, summing all CLEO-c exclusive semileptonic branching fractions gives 
= (6.1±0.2(stat)±0.2(syst))%andE'B(i:'+,i) = (15.1 ±0.5(stat) ±0.5(syst))%. 
These are consistent with the world average inclusive semileptonic branching fractions: 
B{D'^ ^ e+X) = (6.9 ± 0.3)% and B{D+ ^ e+X) = (17.2 ± 1.9)% [13, excluding the 
possibility of additional semileptonic modes of the and with large branching frac- 
tions. 

In summary, we have presented the most precise measurements to date of the absolute 
branching fractions of decays to K^e'^i'e, 7T^e~^i>e, K*^e^Vei and p^e'^Ug, and the first 
observation and absolute branching fraction measurement of ue^Ve- We have com- 

bined these with measurements in Ref. which use the same data and analysis technique, 
to demonstrate that charm exclusive semileptonic decays are consistent with isospin invari- 
ance and to test other theoretical predictions. A comparison of the world average inclusive 
semileptonic branching fractions to the sum of the semileptonic branching fractions in this 
work excludes the possibility of additional semileptonic modes with large branching frac- 
tions. 

The precision achieved in this analysis is consistent with the expected performance of 
CLEO-c. CESR is currently running to collect a much larger ^(3770) data sample. It 
is expected that this sample will result in greatly improved measurements of and 
semileptonic branching fractions, measurements of the decay form factors, which are strin- 
gent tests of LQCD, and the CKM matrix elements Vcs and Vcd 01- 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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